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The effective potential of composite diquark fields responsible for color symmetry breaking in cold 
very dense QCD, in which long-range interactions dominate, is derived. The spectrum of excitations 
and the universality class of this dynamics are described. 
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Recently, there has been considerable interest in the study of the color superconducting phase of cold dense QCD 
ffl-|lOl (for recent reviews, see Ref. fll]]). The color superconducting quark matter may exist in the interior of neutron 
stars, with baryon number densities exceeding a few times the normal nuclear density no — 0.17 fm -3 . Also, such 
matter could be created in accelerators by heavy ion collisions. 

The Ginzburg-Landau (GL) effective action method has been extremely successful in studying ordinary supercon- 
ductivity of metals (l2) . Recently, a similar approach has been utilized in the study of color superconductivity ■ 
However, there the effective action was postulated based on symmetry and renormalization group arguments, and not 
derived from the microscopic theory, QCD. 

Following the original approach of Gorkov jl2| , it would be of a great interest to derive the effective action in color 
superconductivity directly from QCD. In this letter, we make a step in realizing this program and derive the effective 
potential for the order parameter of color superconductivity in cold dense QCD at such high baryon densities when the 
fermion pairing in the diquark channel dominates over that in the chiral one and when long-range interactions 
dominate [0J3]. For this purpose, we will utilize the method of Ref. jn|, which was originally used for the derivation 
of the effective action in quenched strong-coupling QED4 (see also Ref. [Q) and then was successfully applied to 
QED 3 |L5| , quenched QED 4 in a magnetic field |l6| , and to some other models jL7| . 

The crucial feature in the dynamics of cold dense QCD, pointed recently in Refs. |QJ5| (see also Ref. ||), is the 
presence of the long-range interactions mediated by the unscreened gluon modes of the magnetic type. This point 
essentially distinguishes the dynamics of color superconductivity from that in the BCS theory of superconductivity 
in metals. In particular, this makes the derivation of the effective action in color superconductivity more complicated 
than the derivation of the GL effective action from the BCS theory. 

Our derivation of the effective potential in dense QCD will be based on the recent analysis of color superconductivity 
in the framework of the Schwinger-Dyson (SD) equations J?] |J. In this way, we will describe the universality class of 
the dynamics in cold dense QCD and, in particular, get insight into the character of the spectrum of excitations. 

As we shall see below, the universality class of the system at hand is that connected with long-range non-isotropic 
forces, producing a bifermion condensate. We will see that this class resembles (although does not quite coincide 
with) that of quenched QED4 in a constant magnetic field |l8|,[l9| . The scaling law for the order parameter X in these 
two models has the following form: 

X=A eff f(z), /(*)~exp(~), (1) 

and C is some constant. Here z is a generic notation for parameters of a theory, such as a coupling constant, temper- 
ature, the number of fermion flavors, etc.. In QED4 in a magnetic field B, the effective cutoff A. e ff is proportional to 
leJSI 1 / 2 and z is the QED running coupling at the scale ^B] 1 / 2 . In cold dense QCD, which is of main interest here, 
h-eff is proportional to the chemical potential /j, and z is the running QCD coupling constant a s at the scale /z. 

The critical value z c is zero both in cold dense QCD and in QED4 in a magnetic field. This is because in these 
two models, strong interactions are provided by the effective dimensional reduction 3 + 1 — ^ 1 + 1 in the dynamics of 
fermion pairing |]l2| , |l8| . 

One should expect that the long-range interaction in dense QCD leads to the existence of an infinite number of 
resonances in different channels. In particular, as we will see, there is indeed an infinite number of resonances in the 
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channel with the quantum numbers of the Nambu-Goldstone (NG) bosons. It will be shown that this in turn leads 
to a rather unconventional form of the effective potential: it is a multibranched function of the bifermion condensate 
and has an infinite number of local minima. It reduces to the conventional Coleman- Weinberg potential [^0| only in 
the vicinity of the global minimum. 

So, let us consider dense QCD with two light flavors in the chiral limit. The Lagrangian density reads 

Cqcd = i> {i-fDp + M7°) ip ~ \Tr (F^Fn + C gf + C FP , (2) 

where C g f and Cfp are the gauge fixing and the Faddeev-Popov ghost terms. The covariant derivative is defined in 
a usual way, — — ig s A^T A , and p. is the chemical potential. 

Below, it will be convenient to work with the eight component Majorana spinors, $ = ^ ^ c ^ where ip c = C^p T 

and C is a charge conjugation matrix, defined by C~ 1 7 M C = — 7J and C = — C T . The fermion part of the Lagrangian 
density (||) could be rewritten as follows: 

C„, = * ( i7 ' Y " T ° jy , 6 ;_ M/ ) », i, = % + i«,Ai (T-f . (3) 

In order to derive the one-particle irreducible (1PI) effective action for the local composite field 4> c (x) — 

l/2s l: i £ a bc4>a (x)j 5 C (here a, b, c and i, j are the color and flavor indices, respectively), whose vacuum 

expectation value defines the order parameter in the theory, we need to consider the corresponding generating func- 
tional, 

iW{ J c ) =lnj d^d^dA^ exp i J d 4 x (^Cqcd + \ J c e^e abc ^ al 5 C (^) T + c.c. ) . (4) 

When the functional W(J C ) is known, the calculation of the effective action (potential) of interest reduces to performing 
the Legendre transform with respect to the external source J c , 

T(0 C ) = W( J C ) ~ f d 4 x [J c {x)4> c {x) + c.c] , (5) 



where (f> c (x) = (0|<^ c (x)|0) j, and the subscript J implies that 4> c (x) is related to the theory with a source J c . It is 
assumed that the source J c in Eq. (||) is the function of the field (f> c , obtained by inverting the expression, 

<t>o{x). (6) 



SJ c (x) 



For the purposes of calculating the effective potential of the field 4>c, it is sufficient to restrict ourselves to the 
case of a constant (in space-time) external source, J c {x) = Const. In addition, using the freedom of global color 
transformations, it is always possible to fix the orientation of the source in the color space along the third direction, 
i.e., Jx = J% = and J 3 = j 7^ 0. Finally, the baryon symmetry allows us to choose j to be real. 
In the theory with the external source, the inverse of the bare fermion propagator reads 

-1 _ „• ( P + ^7° J 



Go 1 = ~i [ p _ M7 o j, JS = J^-wT 6 . (7) 

Upon neglecting the wave function renormalizations the inverse of the full fermion propagator, G~ 1 , would be 

the same as that in Eq. (^) but with J^ b replaced by ^l b {p) = A(p)£ y £ afc3 7 5 . By inverting it, we obtain the following 
expression for the fermion propagator: 

where 

Ri(p) = (p + M7°) - S (p - 7°St 7 o, (9) 

R 2 (p) = (p - M7°) - 7°^ t 7° (P + ^l )' 1 S. (10) 
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As is clear from the definition of the fermion propagator, A(p) is directly related to the value of the gap in the 
fermion spectrum in the color superconducting phase. At the same time, it is also related to the vacuum expectation 
value of the diquark field. Indeed, by making use of its definition, we obtain 



£ J £ab3 tr 



(Gl2U7 "~ 8 ^ ( 27 r)V§-(b1-M) 2 -A 2 - 



(Note that this expression, up to the change of notations, A — ► <fj7_ = — would remain the same if the gap ansatz 
of Refs. is used. In notation of Ref. |g, A — > Ai.) Therefore, if the solution for the full fermion propagator in 

the problem with an external source is known and the function A(p) is presented, from Eq. ( pd| ) we could also obtain 
the dependence of the diquark field on the source. And, then, it is straightforward to calculate the generating 
functional by integrating the expression in Eq. (^) , 



W{j) r Ao(j) ... ,dj(A 



where, by definition, Ao = A(p)| p= o. 

The gap equation was presented in Refs. There it was also shown that the Meissner effect is of no impor- 

tance for this equation. The further modification of the gap equation for the case of a nonzero external source is 
straightforward, 



1 dq A A(q 4 ) ^ A_ + 2q, f A dg 4 A( g4 ) ^ A (i3) 



Vll + A o P 4 9n Jp* Vll + A o ^4 



where A = (4Tr) 3 / 2 fi/a 5 s ^ 2 ||,^| and a s is the QCD running coupling related to the scale of order fi. This equation, as 
is easy to check, is equivalent to the differential equation, 



P4 A> 4 ) + A> 4) + ^^4==0, (14) 
4 y/pl + A% V 9tt 

along with the infrared (IR) and ultraviolet (UV) boundary conditions, p 4 A'(p 4 )| p4=0 = and A(A) = j, respectively. 
Notice that the dependence of the solution on the source appears only through the UV boundary condition. The value 
of the source itself could be interpreted as the bare Majorana mass. As in Ref. 0, we solve the differential equation 
analytically in two regions p 4 <C Ao and p 4 3> Ao and, then, match the solutions at p 4 = Ao. 
In the region p 4 <§C Ao, the solution that satisfies the IR boundary condition reads 

A(p 4 ) = A J (vJ^) > (15) 

where J n (x) is the Bessel function. In the other region, p 4 3> Ao, the solution, consistent with the UV boundary 
condition, is 

A(p 4 ) = B sin [ V - In — ) + j cos \ V - In — ) . (16) 
V2 p 4 J \2 p A J 

Now, by matching the solutions and their derivatives at the point p 4 = Ao, we get two relations, 

j = A J (z/)cosf|m Aj - A Ji(u) sin In A^J , (17) 
B = A J (u)sm(^ln^j + A J x {v) cos In A^ . (18) 

The first of them relates the value of the gap in the fermion spectrum and the strength of the external source, while 
the other defines the integration constant B in Eq. (|l6|). Now, we can proceed with the calculation of the generating 
functional. By using the relation (O), we calculate the vacuum expectation value of the diquark field, 



V f A dP^iPd _ 16m! aa , (A) = ^ (i 9) 



2 /„ . An2 _i_ A 2 I/ 2 7T 2 ^7T 2 



^ 2 Jo VpI + A 
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with B given in Eq. (||). Then, from Eq. (£|) we obtain the generating functional, 



(20) 



where i? should be considered as a function of j, defined by Eqs. (|17]) and (|18|). After performing the Legendre 
transform, Eq. (EOI) leads to the effective potential V(4>) in the following parametric representation: 



V{A ) 



<HA ) 



» 2 A 2 



2(j »- JlH sin (^ ln ^) 



A 



+4J (z/)Ji(V)cos [vhxj^] -v{j%(y) + Jf(v)) 



8/-t 2 A 



•Ai(")sin ( ^ ln ^ 



+ Ji (i/) cos 



(21a) 
(21b) 



Let us study the properties of this effective potential. In order to determine the vacuum expectation value of the 
diquark condensate, represented by the composite field <j>, we need to know the extrema of the potential in Eq. (I2T 
Thus, we come to the equation dV/dcf) = j(Ao) = 0. By solving it, we obtain an infinite set of solutions for Aq, 



A { n) = Aexp 



Cretan 



2itn 
v 



Aeexp 



3tt 3 / 2 (1 + 2n) 



2 3/2 



n = 0,1,2,. 



which correspond to the following vacuum expectation values of the diquark field: 



%) = (-1)'' 



(22) 



(23) 



Since d 2 V/d4> 2 \j,(i) — (vtt/A^) 2 , we conclude that all the extrema are, in fact, minima. 

It is natural to expect that the potential should also have maxima between those minima. The situation is however 
more subtle: while, as a function of the parameter A , the potential does have maxima, it does not have them as a 
function of 0. Let us describe this in more detail. The first derivative of V with respect to Ao is zero at the following 
maximum points: 



An) 

Omax 



Aexp 



Cretan (#4 



2 /2\ 27T7J 

— arctan — 

v \ v I v 



Aexp 



2nn^ 



1,2,. 



(24) 



However, as is eas y to check, the derivative of the potential with respect to the field <f> at the corresponding 0-points, 
defined from Eq. (21t), is nonzero: it is because the derivative of <f> with respect to A equals zero there. As one 
can see from Fig. |l|, this property is intimately connected with the fact that the potential V{4>) is a multibranched 
(multivalued) function of <f>, and these "maximum" points are sharp turning points at which different branches of the 
effective potential merge. 

As evident from Fig. [j], the physical branch, at which the potential takes the minimal value for a given value of 
(f>, is the first branch, at which the global minimum cf> = (j>^ lies^]- It is interesting that the potential is convex at 
this branch (we recall that the property of the convexity of an effective potential follows from general principles of 
quantum field theory |23|). Moreover, as is seen from Fig. [l], the potential has a fractal structure: after enlargement, 
the higher ("small") branches resemble the first ("large") branch. The whole n-th branch shrinks into the limiting 
point <fi = as n goes to infinity. 

As we will show below, this multivaluedness of the potential is intimately connected with the long-range nature of 
the interaction in the model and implies the existence of many different resonances with the same quantum numbers 
as the NG bosons. 

The global minimum appears at (f>^ . In the vicinity of this minimum the approximate form of the effective potential 
is given by 



1 Another example of a multibranched potential is connected with the (?-term in QCD: the QCD effective potential is a 
multibranced function of the parameter 8 p^ |. The physical branch is again defined as that with the minimal value of the 
potential for a given value of 9. 
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(25) 



i.e., in the vicinity of the minimum, it has the form of the Coleman- Weinberg potential J20|. The region of validity 
of this approximation is given by inequality v\n.{(j>/ (f>^ ) <C 1, and, therefore, Eq. fl25| ) is a very good approximation 
for the potential of the composite field <f> when the coupling is weak or when the value of the field is close to the 
minimum. 

Now, let us discuss how the infinite number of minima in the effective potential (^l| ) determine the form of the 
spectrum of the resonances in the channel with the quantum numbers of NG bosons. It is well known (see, for 
example, Ref. [ p4| ) that, because of the Ward identities for chiral currents, the SD equation for the dynamical fermion 
gap coincides with the Bethe-Salpeter (BS) equation for corresponding (gapless) NG bosons, which are quark-quark 

bound states in the present model. The infinite number of solutions Aq™' (^) for the gap implies that there are 
massless states (which would become the NG bosons) in each of the vacua corresponding to different values of n. The 
genuine, stable, vacuum is that with n = 0. What is the fate of the quark-quark bound states which would be the 
NG bosons in the false vacua, with n — 1,2,...? We will argue below that they become massive, unstable, particles 
there. 

In the chiral limit, there are two free parameters in cold dense QCD: Aqcd and the chemical potential fi, or, 
equivalently, the coupling constant a s (n) and /i. Let us consider the NG composites in a false vacuum, with n=rv- ' > 1. 
In that vacuum, they are massless bound states of fermions with the Majorana mass (gap) being equal to A^°\ The 

transition to the genuine vacuum, with n = 0, corresponds to increasing the fermion gap, Aq" ' 1 — > Aq°\ without 
changing the dynamics: the coupling constant a s (n) and the chemical potential /i remain of course the same. As a 
result of the increase of the mass of their constituents, the square of the mass of these bound states will also increase. 
Therefore they become massive (apparently, unstable) composites in the genuine vacuum^]. 

Thus we conclude that the global minimum Aq of the effective potential indeed defines the dynamical gap (Majo- 
rana mass) of fermions, and all other minima Aq™\ n = 1,2, . . ., manifest the existence of massive radial excitations of 
NG bosons. Notice that, because of the Higgs effect, the NG bosons are "eaten" by the five gluons, corresponding to 
the SU(3) C — > SU(2) C breakdown. All the massive excitations, though, will not be affected by the Higgs mechanism. 

In order to determine the spectrum of these massive excitations, one needs to study the BS equations for massive 
bifermion bound states in dense QCD. This problem is beyond the scope of this letter. However, it is not difficult to 
estimate their masses: since the fermion gap Ag '' is essentially the only relevant dimensional parameter in the pairing 
dynamics, the masses of these resonances should be of the order of the fermion gap. The resonances are unstable, 
although, they might be rather narrow at high density because the coupling constant is weak. The presence of such 
resonances would be a very clear signature of long-range forces in dense QCD. 

Now we come to the description of the universality class of the dynamics in cold dense QCD. The scaling law for 
the order parameter is described by expression ([!]) with X = A ', A e ff ~ [i, and z — a s . The essential singularity 
at a s = is provided by long-range forces. Let us discuss the character of these forces in more detail. 

The gap equation (|l3|) in the absence of the external source can be rewritten in a different form (see Ref. 0), 



4a, r°° dq r A dg 4 A(q 4 ) A_ 

9 J.^ 2tt J_ a 2ir q 2 + q 2 4 + A 2 * \q 4 - p 4 \ ' 



AW^/ ^ r^ ln^-^, (26) 



where the new integration parameter q is the spatial momentum shifted by the chemical potential, q — \q\ — /i. Then, 
it is easy to show that Eq. (|2q) is equivalent to the following Schrodinger equation: 



d?__d?_ 
dr 2 dx 2 



+ A' + U(r,x) )V(t,x)=0, (27) 



where 



2 The quenched strong-coupling QED yields an example of a simpler model with an effective potential having the form similar 
to that in Eq. (J2TJ) ^3p4| . The study of the BS equations in that model shows that there is indeed an infinite number of 
resonances in the channel with the quantum numbers of the NG bosons . Their masses are nearly equal and are of the order 
of the fermion dynamical mass. 
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A(p 4 ) 

_ f) ip4T—ipx (28) 



2a s <., N /" A , , A . 2a s 



C/(r, z) = —^8(x) / dp 4 In = -— ±- [it + 2 si(|Ar|)] <5(z). (29) 

9?r y_ A |p 4 | 9tt|t| 

Here si(z) = — dt sm(t) / 1 is the sine integral function. So, we see that the problem reduces to the Schrodinger 
equation ( p7| ) with a non-isotropic interaction potential presented in Eq. (|2^). This interaction is short range in the 
spatial direction, x, and long-range in the (imaginary) time direction, r [notice that si(z) ~ — cos(z)/z as z — > +oo]. 
It is the latter long-range portion of the interaction that is responsible for the particular scaling law of the order 
parameter as in Eq. 

In some respects, the dynamics in cold dense QCD is similar to the dynamics in quenched QED4 in a constant 
magnetic field |L8[ 19). Indeed, in both these models the dimensional reduction 3 + 1— ► 1 + 1 in the dynamics of 
fermion pairing takes place. This feature and long-range interactions lead to the same scaling law ([!]) for the order 
parameter, which is qualitatively different from that of the BCS type. The form of the effective potentials in these 
two models is also similar (compare Eq. (|2l]) with the expression for the potential in Ref. |lq| ). At the same time, the 
universality class of the system at hand is somewhat different from that in QED4 in a magnetic field: in that model, 
the dynamics is provided by relativistic Coulomb-like forces. As is clear, the difference appears due to the explicit 
breakdown of Lorentz boost transformations in dense QCD by a nonzero chemical potential [notice that there is the 
(l+l)-dimensional Lorentz symmetry in QED4 in a magnetic field]. 

In conclusion, in this letter, we have taken the first step in deriving the effective action in color superconductivity 
of the dense quark matter directly from QCD. In particular, we have derived the 1PI effective potential for the order 
parameter responsible for color symmetry breaking. In this derivation, we used the common assumption that the 
baryon density is high enough, so that the fermion pairing in the color antitriplet channel dominates over that in the 
chiral one. 

The crucial feature in the dynamics of cold dense QCD is the long-range interactions mediated by the unscreened 
magnetic gluon modes Because of these long-range interactions, we argue that the system belongs to a uni- 

versality class that is close to (but not quite the same as) that of quenched QED4 in an external magnetic field 
@@. 

We also argue that the spectrum of the diquark resonances with the same quantum numbers as those of NG bosons 
consists of a very large (infinite in our approximation) number of states with masses of order of the fermion gap. Even 
though these resonances are unstable, they might be relatively narrow at sufficiently high density of quark matter. 
We believe that the presence of such resonances would be a clear signature of the unscreened long-range forces in 
dense QCD. It would be worth studying in detail the properties of these resonances under the conditions produced 
in heavy ion collisions: if the color superconducting phase is ever going to be produced in heavy ion collisions, the 
detection of these diquark resonances might be a crucial piece of information for determining the nature of the phase. 
In further studies, it would be interesting to clarify the properties of resonances in other channels as well as to derive 
the effective potential in the case of intermediate densities, when the chiral and the diquark condensates compete 
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V.E.V. of the diquark field 



FIG. 1. An illustration of the multibranched structure of the effective potential as a function of <j>. The enlargement of the 
higher (second and third) branches is also shown 
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